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A class of algorithms for the Landau gauge fixing is proposed, which makes the steepest ascent (SA) method be 
more efficient by concepts of genetic algorithm. Main concern is how to incorporate random gauge transformation 
(RGT) to gain higher achievement of the minimal Landau gauge fixing, and to keep lower time consumption. One 
of these algorithms uses the block RGT, and another uses RGT controlled by local fitness density, and the last 
uses RGT determined by Ising Monte Carlo process. We tested these algorithms on SU(2) lattice gauge theory 
in 4 dimension with small /3s, 2.0, 1.75 and 1.5, and report improvements in hit rate and/or in time consumption, 
compared to other methods. 
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1. INTRODUCTION 

Gauge fixing degeneracies (existence of Gri- 
bov's copies) are generic phenomena in non- 
abelian continuum gauge theories , and thus there 
exist fundamental problems, e.g., what is the cor- 
rect gauge fixed measure in the path integral for- 
malism. 

Although the foudation of lattice gauge theo- 
ries does not necessitate gauge fixings, they are, 
however, often requeired for field theoretic tran- 
scriptions of their nonperturbative dynammics. 

The principle of the Landau gauge fixing al- 
gorithm is given as an optimization problem of 
some functions along the gauge orbit of link vari- 
ables||. There are many extrema of the opti- 
mization function in general, and all these ex- 
trema points on the gauge orbit correspond to 
the Landau gauge, Gribov copies. The abso- 
lute maximum corresponds to the minimal Lan- 
dau gauge [Q. In order to fix the Landau gauge 
uniquely, the minimal Landau gauge is the most 
favourable goal. 

Local search algorithms were developed by a 
chain of gauge transformations!^,^]. . Since there 
are many Gribov copies, the gauge orbit paths are 
easily captured by these extrema. We call these 
method as the steepest ascent method (SA). In 
such a situation that SA fails to attain an absolute 
extremum, one can try a simple random search[Q 
on the gauge orbit with succeeding SA. 



There is only unique method of Hetrick-de For- 
crandQ (HdeF) aiming at the minimal Landau 
gauge fixing, which is systematic in the sense 
that random trial is not involved. However it 
works successfully only for large (3 samples, rather 
smooth configurations. Thus finding efficient al- 
gorithm for the minimal Landau gauge fixing is 
still an open problem. 

We report results of an attempt in some GA 

type methods in comparison with other methods 

as the simple RGT method. We work on SU(2) 

gauge theory of 8 4 lattice, and define gauge fields 
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-(U x ,n ~ Ul „), then the optimization 



function, fitness, is given as 
1. 



Fu(G) 
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where = GlU Xtfl G x+fl . The fitness F V (G) 
can be viewed as negative of energy of SU(2) spin 
system G sitting on sites, and thus the problem 
is equivalent to finding the lowest energy state 
under the randomized interaction U. 

Straightforward GA strategy was applied to the 
minimal Landau gauge fixing^). Their prelimi- 
nary results of the Landau gauge fixing show that 
the straightforward application of GA is not so 
good as to become a practical method. 

We tested three types of GA methods, which 
are different from each other in how RGT are in- 
corporated in the algorithms. We compare them 
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Table 1 



Search parameters and strategies on SU(2)with (3 — 2.0 performed on DEC a 2100 4/275. 



Method 




Hitrate 


executing time [sec] 


BRGT 


IS on N div = 2 


46/50 


10.9 




IS off N div = 2 


42/50 


11.6 


IRGT 


IS off 


32/50 


11.2 


LFRGT 


Rl = 0.5,i?2 = 0.85 


47/50 


11.6 




Rl = 0.35,i?2 = 0.85 


46/50 


14.1 




Rl = 0.45,i?2 = 0.9 


45/50 


10.3 




Rl = 0.5,i?2 = 0.85 IS on 


44/50 


12.7 


HdeF 




0/50 


7.6 


TRGT 




47/50 


14.7 



with the simple RGT method in performance. 

2. Algorithms 

Our aim is to develop algorithms efficient for 
small /3s. Basic building block for the algorithms 
is RGT. The simple RGT (TRGT) method which 
applies RGT on the whole lattice takes time, 
while it is difficult for SA paths to escape from 
a Gribov copy if the RGT is restricted in blocks 
of too small area on the lattice. 

Thus main features of our algorithms consist in 
how to determine blocks to which RGT is applied. 
Given a Gribov copy link configuration U, the 
following three types of definition of blocking for 
RGT are devised. 

1. The whole lattice is partitioned into Ndi V d 
chequered blocks, where the number of dimen- 
sion, d = 4. Then RGT is applied on white 
blocks and a constant RGT on each black one, 
and vice versa. We call this method as blocked 
RGT (BRGT) method. 

2. We set two parameters Rl and R2, and sites 
i where RGT is applied are chosen according to 
local fitness density, f(i), by f(i) < Rl or R2 < 
f(i). We call this method as local fitness density 
RGT (LFRGT). 

3. Ising spin interaction on randomly chosen 
coarse lattice is defined from gauge spin inter- 
action given by U such that at least one anti- 
ferro interaction should be involved. Through 
Monte Carlo simulation of this Ising system, one 
obtains up-spin blocks B + and down-spin blocks 

Blocks for RGT is given by one of these 



blocks. (3i S i ng is so chosen that size of both blocks 
B± becomes comparable. We call this method as 
Ising RGT (IRGT). 

We use the local exact algoritm|||| f° r SA 
method. Given a Gribov copy U, the SA method 
following RGT in use of one of three blockings, 
brings U to the extrema of fitness by steps of 
gauge transformations. This new copy in the Lan- 
dau gauge is put as an initial copy for the next 
iteration. This sub-procedure is repeated Mu r 
times. The original Gribov copy and the maxi- 
mum fitness copy among Mn r new gauge copies 
are compared in fitness value. If the fitness of 
the new one is higher, the sub-procedure is to be 
started again with this new one as an initial copy. 
Otherwise the process stops, and the initial copy 
is considered as the expected configuration with 
the maximum fitness value. 

In addtion to the above basic algorithm, two 
types of modification are devised as follows: 

1. As the initial copy in the sub-process, the bet- 
ter fitness copy between the current copy and the 
preceding, is always chosen. We call this as the 
inter-selection-on (IS) scheme. 

2. As the initial copy in the sub-process, the 
product of crossing is adopted, where chequered 
block crossing is done between the best fitness 
copy and the second best one among the obtained 
copies so far. We call this as the crossing-on (C) 
scheme. 

We tested these algorithms on SU(2) 8 4 lattice 
with P = 2.0, 1.75, 1.5 and tuned some parame- 
ters, as Ndiv j Rising and Mi tr , with or without IS 
and/ or C scheme. 
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3. Results and Performance 



Our three GA type methods were executed on 
the same set of randomly produced 50 copies 
from a suitably chosen copy from samples, (3 = 
2.0, 1.75, 1.5. The TRGT method and the method 
of HdeF were also tested on the same set. Perfor- 
mance of these methods, hit rate of the minimal 
Landau gauge and average time consumption for 
the gauge fixing, are compared. 

We fix Mitr — 5, and in tests of parameters 
search for f3 = 2.0, we found that BRGT with 
Ndiv = 2 shows a sufficient global search power, 
while BRGT with Ndiv = 4 shows a lower hit 
rate. The IS scheme could be viewed as a kind of 
elitism, and it is known that elitism is a suitable 
strategy when the global search power is avail- 
able. BRGT with IS scheme, Ndiv = 2, shows a 
powerful and efficient search, while IRGT with IS, 
however, does not. For LFRGT, the cut parame- 
ters, Rl and R2, affect its search power. Without 
R2 or with too large R2, even high Rl does not 
work well, nor low R2 without Rl. Since LFRGT 
with both paramters suitabley chosen, achieves 
the efficient search, IS scheme helps. From Ta- 
ble |, our algorithms, BRGT with IS, N dlv = 2, 
and LFRGT with IS, Rl = 0.5 and R2 = 0.85, 
have high hit rates and exhibit good performance 
comparable with TRGT. The HdeF method with 
small (3 is known that it does not tend to the 
maximun fitness [0 . 

The hit rate performance for various methods 
is given in Figure |l|, and the average time con- 
sumption is shown in Table [l| 
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Figure 1. Gribov copies obtained by the vari- 
ous algorithms and the unique copy [ obtained by 
HdeF algorithm. 
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